Thomas-Fermi approach to resonant tunneling in 5-doped diodes 
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We study resonant tunneling in B-5-doped diodes grown by Si-molecular beam epitaxy. A 
Thomas-Fermi approach is used to obtain the conduction-band modulation. Using a scalar Hamilto- 
nian within the effective-mass approximation we demonstrate that the occurrence of negative differ- 
ential resistance (NDR) only involves conduction-band states, whereas interband tunneling effects 
seem to be negligible. Our theoretical results are in very good agreement with recent experimental 
ON ■ observations of NDR in this type of diodes. 
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I. INTRODUCTION 



Resonant tunneling (RT) through semiconductor 
diodes has recently attracted considerable attention be- 
cause of its applications in ultra-high-speed electronic 
devices. Their most remarkable feature is that their 
I-V characteristics show negative differential resistance 
(NDR). Most of these devices are based in double- 
barrier structures, where RT occurs via quasi-bound 
states within the well region in the same energy band 
(conduction- or valence-band). On the contrary, other 
kind of RT devices present interband transitions as n for 
instance, InAs-GaSb or InAs-AlSb-GaSb structures .El In 
the later case, Kane's parameter is not negligible as com- 
pared with bandgaps and, consequently, more elaborated 
band-structures are required to fully account for RT. 

Recently, Sardela et aia have observed NDR at room 
temperature in B-J-doped diodes grown by Si-molecular 
beam epitaxy. A schematic cross section of the stud- 
ied device in shown in Fig. [l[ Experimental NDR peaks 
are seen to appear at about ±0.2 V. However, it is not 
clear whether RT occurs via intraband or interband 
mechanisms: Indeed, the authors argued that prelim- 
inary calculations based on the static solution of the 
Poisson's equation show that the whole structure only 
supports a very shallow quasi-bound state and, con- 
sequently, conduction-band modulation cannot explains 
the observed NDR. Hence, they claimed that the valence- 
band plays a major role in RT of these diodes and there- 
fore two-band Hamiltoniansel are necessary to understand 
the obtained results. 

In this paper we will show that scalar Hamiltonians 
for the conduction-band envelope functions work well to 
explain the results of Sardela et aln Instead of using self- 
consistent approximations to compute the one-electron 
potential due to 5-doping, we will concern ourselves with 
the nonlinear Thamas-Fermi (TF) formulation as intro- 
duced by Ioratti.cl This approach reproduces accurately 
the various subbands densities in periodically (5-doped 



GaAs layers at moderate and high dopingu and it has 
proved itself useful n describing the electronic structure 
in the presence of applied external fields BEI We thus ob- 
tain results showing the existence of several quasi-bound 
states with large lifetimes localized between the two 5- 
layers that successfully explain the I — V characteristics 
previously obtained by Sardela et al. 



II. MODEL 

We consider an electron of effective mass m* in the 
(5-doped diode (see Fig. [y) in the presence of an electric 
field applied along the growth direction. Nonparabol- 
icity effects arc neglected so that m* will be taken to 
be independent of the electron energy. The electron 
wave function and the total energy of the electron are 
given by ^(r,^) = cxp(ik^ • r±)(f>(x) and -E^kjJ = 
E + (H k^/2rn*), where <p{x) is the envelope function. 
Note that we have taken isotropic bands at the X min- 
ima, and, in fact, the band structure of Si is strongly 
anisotropic (m* /m* ~ 5). Hence, we should consider an 
average value m* that would .-be obtained, for example, 
in a measurement of mobilityfl In the effective-mass ap- 
proximation the envelope function satisfies the following 
Schrodinger equation 



h 2 d 2 



2m* dx 2 



+ Vtf{x) -eFx 



{x)=E4>(x). (1) 



Here F is the applied electric field and Vtf(z) is the so- 
lution of the nonlinear TF equation (see Refs. || and [7| for 
details). As usual in scattering problems, the envelope 
function at the bottom electrode is a superposition of in- 
cident and reflected traveling waves, whereas at the top 
one there is only a transmitted wave. Hence standard 
numerical techniquesO can be used to obtain the trans- 
mission coefficient t(E, V) for a given incident energy E 
and a given applied voltage V — FL, L being the length 
of the whole structure. Due to the high doping of the elec- 
trodes, the electric field is assumed to be nonzero only 
within the structure. In addition, these high screening 
effects imply that Vtf(z) also vanishes at the electrodes. 



III. RESULTS 

We have set an effective mass m* = 0.33 and a dielec- 
tric constant e = 11.7 in Si.Q In the absence of external 
fields, the TF potential presents two peaks corresponding 
to the (5-layers, as shown in Fig. 0. It is worth mentioning 
that the maximum value of the conduction-band modu- 
lation is about 20% of the bandgap in Si. This value 
is of the same order of that obtained self-consistently 
in Sb-<5-doped grown by Si-molecular beam epitaxy with 
similar areal concentration of impurities a Assuming a 
similar profile for the valence-band, we are led to the 



conclusion that interband RT should not play a signif- 
icant role, at least at moderate electric fields. Hence, 
the scalar Hamiltonian ([I]) suffices to describe the elec- 
tronic structure of the diode. In addition, since we are 
close to the conduction-band edge, parabolic subbands 
hold, as we assumed previously. To elucidate whether the 
conduction-band modulation can support narrow quasi- 
bound states, we have numerically evaluated the trans- 
mission coefficient at zero bias t(E, 0) as a function of the 
incident energy and results are plotted in Fig. |3|. From 
this figure it becomes clear that four resonances (quasi- 
bound states) appear below the top of the potential. The 
levels are nearly equally spaced in this case: This is easy 
to understand if we notice that the potential profile be- 
tween the two (5-layers is almost parabolic. Interestingly, 
we thus have demonstrated that <5-doped diodes could be 
used to achieve equally spaced peaks in the collector char- 
acteristics, in analogous way to parabolic wells formerly 
proposed by Capasso and Kiehl.tj 

The tunneling current density at a given temperature 
T for the diode sketched in Fig. |]can be calculated within 
the stationary-state model from 

m = nSek^T /- r(E) v)N{E ^ v) dE ^ (2a) 

2-k z Ti Jo 

where N(E, V) accounts for the occupation of states to 
both sides of the device, according to the Fermi distribu- 
tion function, and it is given by 

» R() .JJi*JM) (2b) 

\ 1 + cxp[(E F - E - eV)/k B T\ J 

We show in Fig. [| the computed j — V characteristics 
at T = 77 K and at room temperature. The curves have 
been obtained taking the Fermi energy at the conduction- 
band edge away from the <5-layers. Other values of the 
Fermi energy simply modify the scale of j, keeping its 
shape unchanged. The main NDR feature, also ob- 
servable at room temperature, appears at about 0.22 V, 
which is very close to experimental results of Sardela et 
al.a Note that our computation predicts two separated 
peaks in the j — V characteristics at around 0.22 V, while 
experiments show only a single, broader peak. Inelas- 
tic scattering mechanisms, not included in our analysis, 
are known to cause a broadening p£ the intrinsic level 
width, and the amount of inelastidliJ broadening makes 
the two separated resonances to merge into a single one, 
according to experimental results. Moreover, there exist 
smaller peaks at about 0.10 V and 0.15 V, arising from the 
lower resonances of the well. In experimentsa such small 
peaks are almost unnoticeable, probably because we are 
overestimating its height due to the absence of inelas- 
tic effects in our model which, of course, tend to disrupt 
coherent tunneling. The theoretical peak-to-valley ratio 
of the main NDR feature is 2.8 at 77 , larger than the 
experimental value 1.4 at the same temperature. The 
reason why the theoretical calculations predicts larger 
peak-to-valley ratios is that sequential tunneling is not 
being considered. 



IV. SUMMARY 

We have proposed a simple model to explain RT 
through (5-doped diodes grown by Si-molecular beam epi- 
taxy. The conduction-band modulation is obtained by 
means of the TF approach and the corresponding elec- 
tronic states are described by a scalar Hamiltonian. Our 
results show that NDR effects arc due to conventional 
resonant tunneling process in vertical transport, whereas 
interband tunneling does not give rise to significant con- 
tributions. The obtained results are in excellent .-agree- 
ment with recent measurements by Sardela et aln This 
success suggests that our approach, being very simple 
and computationally inexpensive, may be very useful in 
dealing with semiconductor nanostructures. 
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FIG. 1. Schematic cross section of a B-5-diode grown by 
Si-molecular beam epitaxy, as studied in Ref. ti. The whole 
structure is embebded between two electrodes consisting of 
highly n-doped Si layers. Each J-layer is 8 A width with an 
acceptor concentration 1.4 x 10 13 cm -2 . 



FIG. 2. Modulation of the conduction-band obtained from 
the Thomas-Fermi approach for the structure shown in Fig. hi 
without applied electric field. 



FIG. 3. Transmission coefficient without applied electric 
field as a function of the incident energy. 



FIG. 4. Computed j — V characteristics of 5-doped diodes 
at (a) T = 77 K and (b) room temperature. 



